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Abstract 

We derive a general expression for the gauge invariant mass (mg) for 
an Abelian gauge field, as induced by vacuum polarization, in 1 + 1 
dimensions. From its relation to the chiral anomaly, we show that iuq 
has to satisfy a certain quantization condition. This quantization can 
be, on the other hand, explicitly verified by using the exact general 
expression for the gauge invariant mass in terms of the fermion prop- 
agator. This result is applied to some explicit examples, exploring the 
possibility of having interesting physical situations where the value of 
niQ departs from its canonical value. We also study the possibility of 
generalizing the results to the 2 + 1 dimensional case at finite tem- 
perature, showing that there are indeed situations where a finite and 
non- vanishing gauge invariant mass is induced. 
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1 Introduction 



Some important physical quantities displaying quantization properties, may 
sometimes be represented by means of momentum space integrals which ex- 
hibit their topologically invariant character. Considerable effort has been 
devoted to find these representations, since they are often very useful to 
prove their quantization, as well as their stability under perturbations. 

This is the case, for instance, of the transverse conductivity o in QED3, 
which can be represented as a momentum space integral given by [I]. [2[] 

a = T^eV 1 ^ J ^Tr [d^Sd^Sd^S] , (1) 

where S(p) is the full fermion propagator. As shown in [|TJ, due to the ultra- 
violet behavior of the fermion propagator, this integral reduces to the Kro- 
necker topological invariant which labels the homotopy classes in U 2 (S 2 ), that 
is, a is quantized. A similar representation for the induced Chern-Simons 
coefficient has been obtained in || by using a cubic lattice regularization 
for the Euclidean fermionic action. In this case, the possible values for the 
Chern-Simons coefficient turn out to be labeled by the winding number char- 
acterizing the mapping between the three-dimensional torus in momentum 
space and the normalized quaternion corresponding to the fermion propaga- 
tor: S(p)/^/de~t(S). 

Topology in momentum space has also been advocated in Q], in order 
to discuss the stability of neutrino masses in the Standard Model and the 
spectrum of excitations in effective two dimensional models such as Helium- 
3. On the other hand, in the context of Yang-Mills theories, one of the main 
open problems is that of understanding the non-perturbative generation of a 
mass gap for the gauge fields as a consequence of the dynamics in the infrared 
regime. This problem is of crucial importance for confinement in QCD, and 
for analyzing finite temperature effects. For instance, in 0, a gap equation, 
based on the introduction of gauge invariant mass terms, has been proposed 
and applied to the Yang-Mills theory in 2 + 1 dimensions. 

The aim of this work is to obtain a useful momentum space representa- 
tion for detecting the existence of an induced gauge invariant mass, mg, for 
an Abelian gauge field. In particular, we shall be able to derive a general 
expression for ma, valid for different space-time dimensions and geometries 
in momentum space. 
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The paper is organized as follows. Section 1 is devoted to a short review 
of the well-known relationship between the induced gauge invariant mass and 
the chiral anomaly in 1 + 1 dimensions. In section 2, a general momentum 
space representation for the induced gauge invariant mass is provided. Sec- 
tion 3 is devoted to the applications of the aforementioned representation to 
the case of square lattice geometry and to the case of finite temperature in 
2 + 1 dimensions. Section |5| presents our conclusions. 

2 Relationship Between the Gauge Invariant 
Mass and the Chiral Anomaly in 1 + 1 di- 
mensions 

By gauge invariant mass, mg, we understand the value of the constant ap- 
pearing in a gauge invariant mass term for an Abelian gauge field. In d + 1 
dimensional Euclidean spacetime, a gauge invariant mass term Lagrangian 
C m is given explicitly by 

C m = l -mlA^ u A u (2) 

where 5^ v = S^ u — is the transverse Kronecker 5 in D = d+ 1 dimensions. 

This mass term for is explicitly gauge invariant, although, of course, at 
the price of introducing a non-locality in the action. That is often the reason 
for not including this term in the classical Lagrangian of a standard local 
quantum field theory. In spite of this fact, terms like this naturally arise when 
evaluating radiative corrections to the effective action. In massless QED in 
1 + 1 dimensions, i.e., the Schwinger model [EJ, it is induced by the one- loop 
vacuum polarization graph. Moreover, for such a model, the one-loop result 
is exact, since it does not suffer from higher order corrections. These results 
can also be shown to be related to the chiral anomaly in 1 + 1 dimensions, 
whence the non-renormalization of wiq is inherited. For this model, the non- 
locality can be made good since, in the bosonization approach, the gauge 
field may be written in terms of two scalar fields: A^ = e^dyO + d^ip. In 
terms of these scalars, the gauge invariant mass term becomes a standard, 
local, mass term for the a field. 

In the 1 + 1 dimensional case, the value of ma for a single fermionic flavor 
has the well-known value m G = e 2 /n. In this article, we are interested in 
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deriving general expressions for mg, as a first step to extending the useful 
1 + 1 dimensional result into several non trivial directions. In order to de- 
fine the problem more precisely, we write vtlq in terms of the exact vacuum 
polarization tensor [] Et^ (A;). The Ward identity for IT^ (k), k^Ii^k) = 0, 
allows us to write: 

fl^k) = fl(k 2 ) S^k) (3) 
where U(k 2 ) is a scalar function. Thus mg may also be obtained as: 

m 2 G = lim fl{k 2 ) . (4) 

This definition implicitly assumes gauge and Lorentz invariance, two condi- 
tions that, except explicit indication on the contrary, shall be maintained in 
everything that follows. 

We shall be first concerned with models that can be described by an 
Euclidean action with the following structure: 

S = [ d D xC £ = C F + C G (5) 



where Cf and Cq denote the fermion and gauge field Lagrangians, respec- 
tively. 

Let us begin with a consideration of the 1 + 1 dimensional case, in the 
simplest non-trivial situation of massless QED{\ + 1). This will amount to 
a re-derivation of known results, although we shall present them here from a 
different perspective. The gauge field Lagrangian is 

Cg = -^F^F^y , F^y = d^Ay — d u A^ , (6) 

while fermionic matter is described by a massless Dirac field ip, with a La- 
grangian: 

C F = Hfi + ieA-)^ , (7) 
where we adopted the conventions: 

(7^) f =7/i 75 = 75 , J^Jv = SV + ie t^l5 

1 Thc tilde denotes the momentum space version of an object, whenever it its convenient 
to distinguish it from its coordinate space version. 
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9fj,u — dfiu ^oi — +1 



(8) 



The exact vacuum polarization tensor is, as usual, defined in terms 
of the connected gauge field two-point function G^ w by: 

{A^A u ) conn = , G~l{k) = k 2 5^ u + U^(k) . (9) 

It is easy to see that there is, indeed, a relation between rriQ and the chiral 
anomaly Q(A), since the former may actually be derived from the latter, at 
least for massless QED in 1 + 1 dimensions. To that end, we define Q (A) by 

9,(4) = ; ( 10 ) 

where ( J^) denotes the quantum average of the axial current = —ieip^f^^ip 
in the presence of an external gauge field A^ (we use the cursive A^ notation 
to distinguish it from the dynamical gauge field A^). Hence, the integrated 
form of the anomaly is 

j d 2 xQ{A) = J dFxd^Jl) = j d 2 xe ilv d^J v ) (11) 

where we used the relation: = e^J v , and J M = etp^^ip is the vector 
current. Then, from our knowledge of the chiral anomaly, we may use the 
fact that Q(A) is linear in A to adopt the linear approximation for (J„) in 
flTTp, namely 

(Mx)) = J d 2 yU^{x-y)A v (y) (12) 

with H^ v {x — y) denoting the coordinate space vacuum polarization function. 
It should be noted that the H^ u (x — y) appearing in ( |I~2"D is exact, since only 
an expansion in the external field has been performed. 

Inserting then ( |I2"D into ( ]T0| ) and Fourier transforming, we find that: 

/ d 2 xQ{A) = limfl(fc 2 ) Mme^ikuAuik) , (13) 

assuming that both limits exist. Taking (§) into account, this amounts to: 

J d 2 xG{A) = m 2 G $(A) , (14) 



5 



where $(A) = f d 2 xe llv d il A v is the total flux of e^ v d^A u in Euclidean space- 
time. Equation (|T4]) yields mg in terms of the anomaly. Of course, the left 
hand side in equation (|1^) being the spacetime integral of the anomaly, will 
have the form of a coefficient times §{A), namely: 

J d 2 xG(A) = £$(A) (15) 

thus the content of ([14]) is that the coefficient £ = e 2 /2tt gives precisely the 
value of ttiq. Now, only configurations with e 2 Q(A) G Z may have a finite 
action, and this implies that £ could have only been e 2 /2tt or an integer 
multiple of this quantity, i.e., 

m% = k^ 7 k = 0,1,2,... (16) 

It must of course be possible to prove this quantization, starting from 
the fermionic side of the problem. Indeed, this is the case for the one-loop 
approximation, since we also know that, to that order, the integral of the 
anomaly is in fact the index of the Dirac operator, namely: 

d 2 xQ{A) = ^-(n + -n_) (17) 

where n± denotes the number of zero modes of positive and negative chiral- 
ity in the given background. In our conventions, it is easy to realize that 
n + - n_ = e 2 $(A). 

In the case in which a current-current interaction is introduced, the cor- 
responding effect can be evaluated by means of the bosonization rules, which 
give 

W) = ^-^A v ■ (18) 
it + g 

where g denotes the current- current coupling constant. This relationship 
would imply, through equation flliD , a renormalized gauge invariant mass, in 
agreement with the result of ref.[EJ. 



Also, when the fermions are massive, the right hand side of equation ( Hp 
has to be supplemented with the additional term —2m^}'y^ coming from the 
explicit chiral symmetry breaking. Accordingly, equation (|T3|) reads 



/ d 2 xd tl [Jl) = [ d 2 x (g{A)-2m{^)) = limfl(fc 2 ) lim e lxv ik IJt A v {k) 

(19) 
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Due to the short ranged behavior of the massive Dirac fields, the left hand 
side identically vanishes, which together with equation (f|) gives a zero gauge 
invariant mass tug- In particular, this implies that the contribution of the 
chiral anomaly cancels against that of the explicit chiral symmetry breaking 
term 0. 

3 Momentum Space Representation for the 
Gauge Invariant Mass 

We try, in what follows, to use the fermionic point of view exclusively, in 
an attempt to derive the previous results therefrom, without assuming that 
a loop expansion has been performed. To that end, let us derive an exact 
expression for rriQ in terms of the fermion propagator. We start from the (D- 
dimensional) general expression for IT^, as given by the Schwinger-Dyson 
equations ||, which amount to the diagram of Figure [I], where the white 



p-k 




Figure 1: The exact vacuum polarization graph. 



blobs are included to mean that the lines are full fermion propagators, while 
the black blob represents the full vertex functions. The external legs are of 
course to be truncated, but we have drawn them for the sake of clarity. The 
analytic expression corresponding to Figure |T| is 

IV(*0 = - e2 / T^tr^^r^p-A;)^-*)] (20) 
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where Sp is the momentum propagator and T u the exact vertex function. 
Noting that IL^ may be expressed as in fl3|), we see that 



m 2 G = 11(0) = (D- l)" 1 lim n^fc) 

k— >0 



(21) 



or 



rn 



-e\D - l)- 1 lim | ^ tr [ 7 „ S F (p) T,(p,p- k) S F (p - k)} 



G 



(22) 



Since the theory is gauge invariant, we may relate lim^o r M (p, p — k) to the 
fermion propagator, by using the exact Ward identity: 



lim TAp,p- k) 



,_d_ 

'dp. 



This identity may be used in (|22|) to obtain 

d D p 



rn 



G 



te 2 (D-l)- 1 J 



(27T) 



D 



tr 







(23) 



l^S F {p)—S F l (p)S F (p) 



--^-W^^* 17 ^ 1 " (24) 

which is an expression we shall consider in detail below. It is still formal, in 
the sense that we have not yet made explicit any regularization method. 

Before exploring its divergences in the general case, we note that, to one 
loop order, they are the well-known divergences in the vacuum polarization 
graph of Figure |2|. 

In (1+1) the expression for niQ suffers from both UV and IR divergences, 
as it is clear from fl2~3|), since, when applied to the free Dirac propagator 
S^P = —ip 1 " 1 , it yields: 

In order to make sense of this expression, we exclude a circle of radius e 
around the origin of the momentum plane to avoid the IR singularity, and 



S 



p-k 




Figure 2: The one- loop vacuum polarization graph. 



also use a Pauli-Villars regulator to tame the UV divergences. This amounts 
to defining 



m 



G 



-2e< 



n 



d2 P ^ (—) i 2c 2 
(2tt)2 dp^p 2 ^ 



d 2 p d 



R 



(2tt) 2 p 2 + A 2 ' 



(26) 



where 1Z is the region illustrated in Figure [| and A denotes the mass of 
the regulator field. The momentum P denotes the radius of the integration 
region, and of course P — > oo since the theory is already regularized. 

By a straightforward application of the two dimensional Gauss' theorem, 
we may convert these integrals of divergences into fluxes of radial vector 
fields. This procedure yields 



2e 2 



m 



G 



-2n £ - 

e 



2tx 



e 2 + A 2 



(27) 



which in the e — > limit becomes: 



m 



G 



e 

7T 



(28) 



as it should be. The procedure we have followed has a simple electrostatic 
analogy: mg is given by the integral of the divergence of an 'electric field', 
so by following that analogy it is proportional to the total 'electric charge'. 
Because of the regulator, though, only the charges at the origin are relevant; 
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Figure 3: The integration region for the integrals in (|26| ). 



namely, the subtraction due to the UV regulator leads to: 




Jc(e) 



(29) 



where the integral is taken along a small curve of radius e enclosing the 
origin, denotes the outer normal to C(e). In what follows we shall argue 
that this kind of expression can be generalized to the full theory. Indeed, we 
can always say that the role of any UV regularization will be to modify the 
large momentum behavior of the propagator, in such a way that all the points 
with infinite momentum can be identified. For the Pauli-Villars case, this 
can be shown to hold simply by combining the contributions of both fermion 
propagators into one integral, to define a regularized propagator. Namely, 



which contains the divergence (in momentum space) of a vector field which 
decreases as ~ p~ A at infinity. 





d 2 p d [ A 2 



(30) 



(2tt) 2 d Vil \_p\p 1 + A 2 ) 
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This will be one of the properties we shall demand of a regularization 
at any finite order of the loop expansion; namely, the points at infinity in 
the momentum space may, from the point of view of this calculation, be 
identified. The other condition is that, in the small momentum region, i.e. 
momenta much smaller than the cutoff, the behavior of the propagator should 
be the same as for the unregularized propagator. Thus only the small region 
around zero may contribute. This, in turn, will produce a non-vanishing 
answer only when the field is massless, as we shall see now. 

Assuming a standard Dirac action for the fermions, the general form of 
the fermion propagator in momentum space is, of course, 

Sp'ip) = tA{p)p'+B{p) (31) 

where A and B are real functions that depend only on the scalar p 2 . Then 
an application of ([24] ) yields: 

2 • 2 f d2 P d . 

m r = —te / -. — r- — — tr 

J (2tt) 2 d Pfl 



7/i 



Mp) I 
A 2 (p)p 2 + B 2 (p) 



W? tr(/) 1, d ' K 



A 2 (p)p 2 + B< 



(32) 



where tr(J) counts the dimension of the Dirac algebra representation. Eval- 
uating the integral along C(e), 



711 



G 



27T 



tr(J) 



A(e)e' 



A 2 \ 



e e 



+ B 2 ( 



(33) 



Since e — > 0, we may conclude from here that, if the fermion is massive, 
B(0) 0, then rrtQ = 0. This is indeed the case for the massive Schwinger 
model. Regarding the massless case, we obtain: 



m 2 r = — tr(J) — ^- = — tr(J) 
G 2ti K J A(0) 2tt V ; 



(34) 



where we have used the normalization condition A(0) = 1, which holds at 
any finite order of the loop expansion. This condition fixes the residue of the 
perturbative electron propagator at the pole p 2 = 0. Strictly speaking, this 
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normalization will change in an interacting theory. Indeed, the fermion prop- 
agator can in general, be rewritten by using the spectral representation ||: 

SHp) = fv ^ ( ^? W (35) 
Jo P + ^ 

where p\ and p 2 are real functions. It is then clear, by the linearity of this 
expression, that there will be a finite gauge invariant mass if there is an 
isolated pole at zero momentum, namely, if we can write: 

s FiP) = * + r ^ zMqi+^i (36 ) 

P Jm? P + A* 

where m? is the multiparticle threshold. In this situation, we would obtain, 

e 2 1 e 2 

m 2 r = — tr(J) — — = — tr(/)Z 2 , (37) 
G 2n K J A(Q) 2tt v ; ' v ; 

where Z 2 is a constant, smaller that 1 because of the spectral condition: 

dp 2 Pl (p 2 ) . (38) 



Equation (|37|) may seem to contradict our remarks on the relation between 
tuq and the anomaly, for example (|i~6|) . The resolution of this apparent 
paradox is that the anomaly, as we understood it in (0), is defined in terms 
of the divergence of the current. The matrix elements of this current, when 
evaluated through the reduction formulae, will require the introduction of 
a Z 2 factor {Z^ for each fermionic field), thus the proper relation that 
generalizes ([16]) is: 

e 2 

m 2 G = k — Z 2 , k = 0,1,2,... (39) 

The original expression for mg required a regularization procedure in 
order to be well defined. However, being determined by the low momentum 
behavior of the propagator, ma should be independent of the regularization. 
For example, for a generalized Pauli-Villars regularization in the one-loop 
case, one might include a set of iV regulator fields, and define: 

N 1 

n=0 r 
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where M n denote the regularizing masses, and the index n = is reserved 
for the unregularized propagator, which has M = and Cq = 1. The 
masses M n and the coefficients C n for n > are chosen in order to verify 
the desired behavior in the UV. However, since only the behavior around 
zero momentum is relevant, and all the regulators are massive, then just the 
unregularized massless propagator contributes. 



4 Applications to Different Geometries and 
Space-Time Dimensions 

4.1 The case of the 1 + 1 dimensional square lattice 

It should be clear from what we have said above that the geometry of the 
momentum space is crucial in the problem of evaluating tug- An interesting 
example of this is the case of a lattice regularized theory. Indeed, assuming 
that the coordinate space has been discretized to an (infinite) square lattice 
with lattice spacing a, the momentum space becomes a (continuous) torus. 
The lattice points are defined then as the set of points = at^ with ^GZ 
for n = 0, 1. 

The naive (i.e., no Wilson term) lattice fermion propagator for a massive 
Dirac fermion then becomes: 

Sf(p) = ■ ttAt (41) 

where C^(ap) = sm(ap^), and the p^ are continuous variables (because the 
lattice is infinite), in the first Brillouin zone, namely: 

-7T < P M < 7T , (1 = 0,1. (42) 

Then, the expression for mg is: 



h ( 27r ) % 



dp d 



Qiap) + am 



(43) 



where B denotes the first Brillouin zone. The fact that we are dealing with a 
momentum torus is obvious because of the periodicity of the functions C^, a 
property which is preserved of course for more general propagators. But the 
functions C M will simultaneously vanish for the (po,Pi) values in the set 

{(0,0), (0,vr), (tt,0), (7r,7r)} , (44) 
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which contain not just the origin but also three other unwanted points. The 
behavior of the lattice propagator is identical for each of these points, and 
equal to the one in the continuum propagator. Indeed, C^ap) ~ a{p — 
where q a is any of the points where both C^'s vanish. 

Then, when m = 0, an application of the Gauss law on the torus, with the 
circles around the four poles of the propagator excluded yields four times the 
contribution of the physical pole at zero. This is because of the contribution 
of the 'charges' at the doublers. This is just a manifestation of the 'doubling' 



problem of lattice fermions, related to the Nielsen-Ninomiya theorem [|TT 
Thus, 



e 2 

m 2 G ]lattice = 4— . (45) 



7T 



4.2 QED% at Finite Temperature 

Another circumstance where the structure of momentum space allows for 
the emergence of a non-trivial gauge invariant mass, is the case of QED 
in 2 + 1 dimensions at finite temperature in the presence of 'large' A Q field 
configurations ||12|| . Parity conserving versions of QED(2 + l), like the model 



introduced by Dorey and Mavromatos |13| have been extensively studied as 
quantum field theory models at finite temperature [)12| . Most of what we 
shall say about this here, however, holds true both for the parity conserving 
and the parity breaking cases. 

We now discuss QED(2 + 1) with regards only to one of its aspects, 
namely, its gauge invariant mass. We recall that in finite temperature quan- 
tum field theory Lorentz invariance is lost. Indeed, in the Matsubara formal- 
ism, which we shall adopt, the time coordinate runs from to j3 — ^. This 
lack of Lorentz invariance implies that one shall, in principle, have different 
masses for the spatial and temporal components of the gauge field. Namely, 
the natural extension of (0) to this case would be 

C m = ^m 2 el A (x,t)A (x,t) + -m 2 mag Aj(x,T)5j- k A k (x,T) (46) 

where 5f k = 5j k — is the transverse Kronecker 5 for the 2 spatial dimen- 
sions, and r is the imaginary time. The two components m e ; and m mag are 
the electric and magnetic masses, which can of course be different for T ^ 0. 
These two masses play the role of components of a 'mass tensor' for A^. We 
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shall deal exclusively with the magnetic mass. It is not difficult to apply a 
similar derivation to the one used for the T = case, to obtain 



m li 



£ Swr^ blSFM] (47) 



where the sum runs over the Matsubara frequencies u n = (2n + l)vrT, and 
Sp is the finite temperature fermion propagator. 

For the free, massless fermion propagator in the presence of a large A , 



we have 12 



S F ( P ,n) = , . 1 (48) 

where u n = u n + eA . A is assumed to be constant, a fact what can always 
be achieved by a small gauge transformation. The value of this constant is 
of course determined by the quantity j drAo(r), where Aq(t) is an arbitrary 
time dependent configuration. Then, 



2 • 2 

m ol = —te 



(49) 



which looks like a series of 1 + 1 dimensional contributions. Of course, when- 
ever u n = — eAo, we have a non- vanishing contribution, of the same kind as 
those appearing in 1 + 1 dimensions. That condition on Aq amounts to 

(1tA (t) = (2n + 1)tt . (50) 
Thus, we may write the result for m 2 el as: 



in. 



e 2 



2irf3 



+oo „p 

tr(/) V S[e / drA (T) - (2n + 1)tt] . (51) 

l=-oo Jo 



The (3 dependence should have been expected by dimensional analysis, 
since in 2+1 dimensions e 2 has the units of a mass. Regarding the existence of 
particular points where the magnetic mass is generated, the physical reason 
for that is that when the condition for Aq is met, dimensional reduction 
occurs because there is a massless mode, and the reduced model is then 
tantamount to a Schwinger model. 
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5 Conclusions 



In this work a useful momentum space representation for the gauge invariant 
mass has been obtained and applied to different situations, namely: massless 
and massive QED in 1 + 1 dimensions, the 1 + 1 dimensional square lattice 
and QED3 at finite temperature. The key ingredient for this representation 
is the validity of gauge invariance expressed by the Ward identity, and of 
course its consistency with the dynamics defined by the Schwinger-Dyson 
equations. 

It should be clear from fl24|) that the value of n? G depends entirely on 
the infrared behavior of the fermion propagator, resulting indeed from the 
infrared dynamics. In particular, the momentum representation for the gauge 
invariant mass in 1 + 1 dimensions turns out to be related to the chiral 
anomaly. This relationship suggests that mc should possess some kind of 
stability against perturbations. For instance, the introduction of a current- 
current interaction amounts to a smooth change |J] of the coefficient A(0) in 
equation fl34|) , while preserving the vanishing of B(p) = 0. This means that, 
apart from a normalization factor, w? G is still finite and non- vanishing. 

A further comment we would like to add, and which is related to the 
previous discussion is one concerning the meaning of formula The 
main difference between our momentum space representation for m G and 
similar representations fl], || |J lies in the presence of the factor A(0). This 
is related to the fact that in our case the formula cannot be written just 
in terms of the SU(2) projection of the fermion propagator as in the other 
cases. Therefore, when including interactions niQ is not protected from being 
renormalized. In this sense, the quantized gauge invariant mass of the 1 + 
1 dimensional Schwinger model is not stable against interactions. In the 
present case, what is stable against perturbations is the ratio between mg 
and Z2. If perturbations are included in the massive Schwinger model or in 
the massless case, wig will however still be tuq = or mc 7^ 0, respectively. 

We finish with a comment and outlook on the possible implications of our 
results to higher dimensional systems at zero temperature, where dimensional 
reduction does not occur. It should be noted that the main difficulty is the 
fact that rriQ, as given by (p4|), seems to be not interesting for the higher 
dimensional case, because of the IR behavior of the fermion propagators. 
One possible way out of this, could be to consider non-standard fields like, 
for instance, dipole fields. This, however, should require a re-derivation of the 
main results, since the structure of the model will, in principle, be different 
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from the standard minimal coupling case. Results on this will be reported 
elsewhere. 
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